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Abstract 

The one-dimensional shock structures of magnetosonic waves (MSWs) propagating in a dissipa- 
tive quantum plasma medium is studied. A quantum magnetohydro dynamic (QMHD) model is 
used to take into account the quantum force term due to Bohm potential and the pressure-like 
spin force term for electrons. The nonlinear evolution (Korteweg de-Vries-Burger ) equation, de- 
rived to describe the dynamics of small amplitude MSWs, where the dissipation is provided by 
the plasma resistivity, is solved numerically to obtain both oscillatory and monotonic shock struc- 
tures. The shock strength decreases with increasing the effects of collective tunneling and increases 
with increasing the effects of spin alignment. The theoretical results could be of importance for 
astrophysical (e.g., magnetars) as well as for ultracold laboratory plasmas (e.g., Rydberg plasmas). 
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The subject of quantum plasmas have received a great attention in investigating vari- 
ous collective quantum effects in plasmas [see e.g. Refs. 1-14]. The collective motion of 
Fermi particles in a magnetic field thus gives rise a natural extension to the classical theory 
of magnetohydrodynamics (MHD) in terms of the well-known quantum magnetoplasmas, 
which have potential applications in astrophysical plasmas, such as pulsar magnetospheres, 
magnetars. Moreover, the motion of particles with spin properties become important in 
strong magnetic fields as a probe of quantum physical phenomena [15-17] in the laboratory 
plasmas. Many of these studies are motivated on a single particle properties. It is thus ex- 
pected that the collective spin effects can influence the propagation characteristics of waves 
in a strongly magnetized quantum plasma [18-20]. Moreover, recent progress in producing 
Rydberg plasmas may give rise to an interesting experimental evidence for the dynamics of 
quantum plasmas. However, in such magnetized plasmas the thermal energy of the particles 
can be very small compared to the typical Zeeman energy of the particles. Recent investiga- 
tions indicate that the spin properties of the electrons and positrons can lead to interesting 
collective effects in quantum magnetoplasmas [19]. More recently, it has been shown that the 
electron spin 1/2 effect significantly modifies the dynamics [21] and modulational instability 
domain [22] of magnetosonic solitary waves and the collective effects in strongly magnetized 
plasmas [23]. 

There has also been much interests in investigating strucures and dynamics of shock 
waves in various quantum plasma media [24-27]. The dynamics of classical shocks is gov- 
erned by a Korteweg de-Vries-Burger (KdVB) equation. A stationary solution of the latter 
can be represented as an oscillatory shock. However, when the dissipation overwhelms the 
dispersion and when the dissipative effect is in balance with the nonlinearity, we indeed have 
the possibility of monotonic shock waves. Unlike the classical fluids, quantum plasmas typ- 
ically exhibit dispersion due to the collective tunneling associated with the Bohm potential 
instead of dissipation. For this reason, even a quantum shock propagating with constant 
velocity in a uniform medium does not exhibit a stationary structure. Transition from initial 
to compressed quantum media occurs in the form of a train of solitons propagating with 
different velocities and with different amplitudes. 

In this letter, we derive a governing equatuion that describes the dynamics of magne- 
tosonic waves (MSWs) in a quantum electron-ion plasma. The governing KdVB equation 
contains both dispersive term due to Bohm potential and the dissipative term due to plasma 
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resistivity (neglecting other effects viz., thermal conduction, viscosity etc.), and also the 
pressure-like spin quantum force. Still, when the normalized zeeman energy~ 1 and the 
plasma resistivity is small, we can recover monotonic transition of the oscillatory shocks. 
The stationary shock solutions exist for the Mach number> 15. The effects of collective 
tunneling and spin alignment influence the strength of the shocks. 

The basic set of equations governing the dynamics of the magnetosonic waves in a quan- 
tum plasma reads [21,22] 

dB d ,„ , d 2 B n 

where B is the magnetic field along the z-axis, i.e. B = B(x,t)z, normalized by its equi- 
librium value B ; p = (m e n e + miUj) = p(x,t) is the total mass density normalized by 
its equilibrium value po; and v = (m e n e v e + mjnjVj)/p = v(x,t)x is the center of mass 
fluid velocity normalized by the Alfven speed Ca = ^/Bq/ p p . The space and time vari- 
ables are normalized by, respectively, Ca/^o, and the ion gyroperiod cu^ 1 = (eBo/rrii) -1 . 
Here n e (rii) is the electron (ion) number density, m e {rrii) is the electron (ion) mass, v e (vj) 
is the electron (ion) fluid velocity and e is the magnitude of the electron charge. Also, 
(3 = 2c 2 {m e /mi)oj'l i \ 2 c /C\ , where = c/u c = h/2m e c is the Compton wavelength, uc 
is the Compton frequency, c is the speed of light in vacuum, h is the Planck's constant 
divided by27r, c s = ^Jk B (T e + T,j)/mi is the sound speed, where T e (Tj) is the electron (ion) 
temperature, and k B is the Boltzmann constant. Moreover, 7 = i]u} ci / p C A , where 77 is the 
resistivity, v\ = kBT e /miC A = ( y l/e)pBBo/m i C A with \x& = eh/2m e is the Bohr magneton 
and e = psBo/ksTg is the temperature normalized Zeeman energy. Note that in Eq.(2) 
[the second term on the right hand side] we have used the isothermal equation of state 
for electrons as P e = kBn e T e for one-dimensional magnetosonic wave propagation across 
Bq. One can also use the equation of state for electrons as P e = m e Vp e nl/3nl for one- 
dimensional propagation [3] or P e = m e Vp e n^ 3 / 5n^ 3 in three-dimension, assuming a local 
zero-temperature Fermi distribution [28] . The last two terms in the right-hand side of Eq. (3) 
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are due to the effects of collective tunneling and spin alignment, respectively, and in Eq.(3) 
we have neglected the inertial term. 

In order to investigate the dynamics of MSWs, we employ the standard reductive pertur- 
bation technique (RPT) with the following stretching 

i = e l /\x-v Q t) J T = e*'% (4) 

where e is a small expansion parameter and vq is the wave phase velocity normalized by Ca- 
The dynamical variables are expanded as 

p = l + ep 1 + e 3/2 p 2 + e 2 p 3 + 

v = ev 1 + e 3/2 v 2 + e 2 v 3 + (5) 

B = l + eB 1 + e 3/2 B 2 + e 2 B 3 + .... 

Now, substituting the expressions [Eqs.(5)] into the Eqs. (l)-(3) and collecting the terms 
in different powers of e we obtain in the lowest order of e 

Pi = B u vi = v (B 1 - 1), (6) 

together with the linear dispersion relation: 



Vq — J 1 + c 2 2" (2 tanh e — e sec h 2 ^) . (7) 

V V B 

From the next order of e, we obtain 

7 dB\ dB 1 
p 2 = B 2 -\ — , v 2 = v B 2 + 7—— - v (8) 



and 



Vq + V V 2 



1 2" (tanh e + e sec h 2 e) 



V B 



B 2 + [ c 2 - 4rtanhe j p 2 (9) 



Inserting Eq.(8) into the Eq.(9) we obtain 



7 ^ 2 -c 2 + ^tanh^^ = (10) 

Since the second factor in Eq.(10) is non-zero by means of Eq.(7) and also dB 1 /d£ ^ 0, 
7 should be at least of the first order of e, so that r ydB 1 /dt; becomes of the order of e 2 , and 
it will be included in the equations for the order of e 2 . Collecting the terms in powers of e 2 
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and eliminating the quantities P3, 1*3 [the coefficient of B 3 becomes zero by Eq.(7)] we obtain 
with the help of Eq.(6) the required KdVB equation 



db db ^b 8 2 b n 

]Tt + p % + q W + r W = °> 

where b = B x and the coefficients P, Q and R are given by 

(12) 



1 

2^o 



3 - vl + 2c 2 s r (8 tanh e + 7e sec h 2 e - 2e 2 tanh e sec h 2 e) 



Q = —7—, R = 7-7 c — M 5- tanhe . (13) 

4V 2w 2 V v 2 B J 

Note that the spin quantum effects are embedded in all of P, Q and R, whereas the dis- 
persion due to quantum diffraction and dissipation due to plasma resistivity are in Q and R 
respectively. We now numerically solve the Eq.(ll) directly in order to obtain nonstation- 
ary shock solutions. In the numerical scheme the KdVB equation (11) is advanced in time 
with a standard fourth-order Runge-Kutta scheme with a time step of 10~ 4 s. The spatial 
derivatives are approximated with centered second-order difference approximations with a 
spatial grid spacing of 0.2 m. The profile of the oscillatory shock solution of Eq.(ll) for the 
parameter values B = 0.14T, T e = 0.09K, e = 1.04, n = 10 30 m- 3 , A = r]/fi = 0.001 is 
shown in Fig.l. The train of oscillations propagates together with the shock with the same 
velocity. As the role of spin force increases, the shock strength decreases and the oscilla- 
tions ahead the shock becomes less in number, in which the first few oscillations are very 
close to the magnetosonic solitons. The oscillations decay quite slow as the role of quantum 
diffraction increases. The plot with A = 0.01 shows the monotonic transition (Fig. 2) from 
the oscillatory shocks shown in Fig.l. Increasing further the role of quantum effects, we can 
not observe the oscillatory shock transition from the Fig. 2. 

The stationary solution of Eq.(ll) can also be obtained by transforming to the moving 
frame of reference ( = £ — Vt = y/e(uj ci /CA) [x — Ca(v + eV)t\ . The KdVB equation then 
reduces to the following system: 

db da If P,, P\ 

TC a 'TC~Q{ Ra+ 2 b - Vb + v -l) < 14 > 

The system of equations (14) has two singular points, namely (1,0) and (2V/P — 1,0) 
which are stable node or focus according as R 2 + AVQ ^ (since Q < 0). A stable focus 
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corresponds to an oscillatory shock (Fig.3) (dispersion dominant), while a stable node gives 
rise monotonic shocks (Fig. 4) (dissipation dominant). The shock strength is given by 

[^ max = e (^-l)=^(M-l)- e , (15) 

where we have defined the shock Mach number M as the ratio of the velocity Ca(vq + eV) 
of the nonlinear magnetosonic wave to the linear wave velocity CaVq by 

M=l + e— . (16) 

Numerical solutions of Eq.(14) are shown in Figs. 3 and 4 for the same parameter values 
as in Figs.l and 2 respectively, but for V — 20. We find that for these sets of parameters 
shock solutions exist for V > 15. As the value of V increases, the number of oscillations ahead 
of the shock increases with decreasing the shock amplitude near ( = 0. Also, as the value 
of the zeeman energy e decreases, the shock strength increases and the oscillations decay 
quite slowly forming long wave train, while for large value of e, the oscillations decay quite 
fast. Numerical simulation also reveals that the shock strength decreases with increasing 
the particle numer density and decreasing the electron temperature, while it increases with 
the strength of the ambient magnetic field. 

To conclude, we have investigated the effects of quantum tunneling and spin alignment 
on the magnetosonic shock structures in a dissipative quantum plasma medium. The numer- 
ical solutions of the KdVB equation exhibit both stationary and nonstationary oscillatory/ 
monotonic shock solutions in the quantum regime. Such significant modifications of the 
shock structures in our quantum plasma are completely a new feature relevant for astro- 
physical and ultracold laboratory plasmas. 
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